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Abstract 

We study the exclusive decays of B — > K*£ + £~ within the framework of the perturbative 
QCD (PQCD). We obtain the form factors for the B — > K* transition in the large recoil 
region, where the PQCD for heavy B meson decays is reliable. We find that our results for 
the form factors at q 2 = are consistent with those from most of the other QCD models in the 
literature. Via the decay chain of B — > K*(K7r)£ + £~ , we obtain many physical observables 
related to the different helicity combinations of B — > K*£ + £~. In particular, we point out 
that the T violating effect suppressed in the standard model can be up to 0(10%) in some 
CP violating models with new physics. 

Key words: B decays, perturbative QCD, CP violation 



1 Introduction 



There has been an enormous progress for flavor physics since the CLEO observation |TJ of 
the radiative b — ► sj decay. Recently, the decay modes of B — > K£ + £~ (£ = e, jj) have been 
observed || at the Belle detector in the KEKB e + e~ storage ring with the branching ratio of 
Br(B K£ + £~) = (0.75igjf ± 0.09) x 10~ 6 , while the standard model (SM) expectation is 
around 0.5 x 10~ 6 [Q. We remark that the decay has not yet been seen at the BaBar detector 
in the PEP-II B factory [|[]. Experimental searches at the B- factories for B — > K*£ + £~ are 
also within the theoretical predicted ranges |J. It is known that the study of flavor change 
neutral currents (FCNCs) in these B decays provides us with information on not only the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix elements |J in the SM but also 
new physics such as supersymmetry (SUSY). 

On the other hand, via B decays such as B — > J/\& K, we can test whether the unique 
phase in the CKM matrix is indeed the origin of CP violation (CPV). In general, CP asym- 
metries (CPAs) in B decays are defined by acp oc V — T and Acp{t) oc T(t) — T(t), called 
direct CPA or CP-odd observable and time dependent CPA, respectively. The former needs 
both weak CP violating and strong phases, while the latter contains not only a non-zero 
CP-odd phase but also the B — B mixing. We note that the present world average for a^p 
is 0.79 ± 0.12 0, § comparing with the SM prediction of 0.70 ± 0.10 [§. In the decays of 
B — K*£ + £~, CPAs such as acp are small even with weak phases being 0(1) due to the 



smallness of strong phases ]T0 



To study CPV in B — > K*£ + £~ , one can also define some T-odd observables by momen- 
tum correlations, such as the well known triple momentum correlations 0. These observables 
do not require strong phases in contrast to the CPA of acp- In the absence of final state 
interactions, these T-odd observables are T violating and thus CP violating by virtue of 
the CPT theorem. In the decays of B — > K*£ + £~ (£ = e,fi, and r), the spin s can be 
the polarized lepton, s#, or the K* meson, e*(A). For the polarized lepton, since the T-odd 
polarization is normally associated with the lepton mass, we expect that this type of T vio- 
lating effects is suppressed and less than 1% for the light lepton modes |ll|]. Although the r 
mode can escape from the suppression, the corresponding branching ratio (BR) of O(10~ 7 ) 
is about one order smaller than those of e and /i modes. In this paper, we concentrate on 
the decay chain of B — > K*£ + £~ — > Ktt£ + £~ and give a systematic study on various possible 
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physical observables, especially the T-odd ones. 

It is known that one of the main theoretical uncertainties in studying exclusive hadron 
decays arises from the calculations of matrix elements. At the large momentum transfer (q 2 ) 
region, Lepage and Brodsky (LB) Jl2[ have developed an approach based on the perturbative 
QCD (PQCD). In the LB formalism, the nonperturbative part is included in the hadron 
wave functions and the transition amplitude is factorized into the convolution of hadron 
wave functions and the hard amplitude of valence quarks. However, with the LB approach, 
it has been pointed out that the perturbative evaluation of the pion form factor suffers a non- 
perturbative enhancement in the end-point region with a momentum fraction x — > [|TB|]. If 
so, the hard amplitude is characterized by a low scale and the expansion in terms of a large 
coupling constant a s is not reliable. Furthermore, more serious end-point (logarithmic) 



singularities are observed in the B — > n transition form factors |L4], [15[ from the twist- 
2 (leading-twist) contribution. The singularities become linear while including the twist-3 



(next-to-leading twist) wave function fl6 |. Because of these singularities, it was claimed that 
even at the low q 2 form factors are dominated by soft dynamics and not calculable in the 
PQCD 0. 

Following the concept of the PQCD, if the spectator quark inside the B meson with a 
momentum of 0(A), where A = Mb — vn^ and mj is the b quark mass, wants to catch up 
the outgoing quark with an energy of 0(Mb/2) to form a hadron, it should obtain a large 
energy from b or the daughter of it. That is, hard gluons actually play an essential role in 
the B meson with large energy released decays. Therefore, relevant decay amplitudes should 
be calculable perturbatively. It is clear that to deal with the problem of singularities is the 
main part of the PQCD. In order to handle these singularities, the strategy of including 
kf, the transverse momentum of the valence-quark [18|], and threshold resummation |19 , pO| 
have been proposed pi| . It has been shown that the singularities do not exist in a self- 
consistent PQCD analysis [EI]. In the literature, the applications of this PQCD approach to 
the processes of B -> PP, such as B -»• Kir pf, B -> irir @, B -> KK [H, B -> Krj^ 



[2§ and B s -> KK [|6|, as well as that of B -> VP, such as B -> </>tt [gj, B — > cf>K 
2"Bfl, B — » p(uj)n [^9| and B — > p(w)K |3(J, have been studied and found that they are 



consistent with the experimental data. In this paper, to calculate the matrix elements of 
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relevant current operators, we adopt the PQCD factorization formalism as 

~d 2 b 



(V\O k \B) = J[dx 

x S t ({x})e 



$y(x 2 , 62) T k ({x}, {&}, M B )<Z> B ( Xl , h 



Air 

-S({x},{b},M B ) 



where $y ($b) is the wave function of V (B) meson, 7\ is the hard scattering amplitude 
dictated by relevant current operators, the exponential factor is the Sudakov factor j32fl , 
and St(x) |3~3], |3~1] expresses the threshold resummation factor. 



The paper is organized as follows. In Sec. II, we study the form factors of the B — > K* 
transition in the framework of the PQCD. In Sec. Ill, we write the angular distributions and 
define the physical observables for the decays of B — > K*£ + £~. In Sec. IV, we present the 
numerical analysis. We also compare our results with those in other QCD models. We give 
our conclusions in Sec. V. 

2 Form factors in B — ► K* 

For decays with the B — > K* transition, the B meson momentum p\ and the K* meson mo- 
mentum p 2 and polarization vector e in the B meson rest frame and the light-cone coordinate 
are taken as 

M B , - v M B . 2 2 - . 

Pi = -^=(l,l,0j_), p 2 = -pf(V ,r K *,a±), 

= (r,\-r 2 K *,0 ± l e T (±) = -Uo, 0, 1, ±i) (2) 

with 1] « 1 — q 2 /M B and tk* = Mk*/Mb, while those for the spectators of B and K* sides 
are expressed by 

h = (o, Xl ^=, ki±j and k 2 = (x2^^rj,0,k 2 i\ , (3) 

respectively. In our calculations, we will neglect the small contributions from m u ^ s and A as 
well as M|-» due to the on-shell condition of the valence-quark preserved. From the results 



in Ref. [35], the K* meson distribution amplitude up to twist-3 can be derived as follows: 



+M K *[I] lj( j ) Ux)h 
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(lf(p,6 r )|sW i d(0) I |0> = ^J^xe^iM^lM^U^ + iAT fa<&.{x) 

+j^e^[^] l ^nZ.<^ K .{x)}, (4) 

where n_ = (0, l,0j_) and <pK*{x) and are the twist-2 wave functions for the longi- 

tudinal and transverse components of the K* polarization, respectively, while the remaining 
wave functions belong to the twist-3 ones with their explicit expressions given below. 

2.1 Power counting 

To show the B — > K* form factors, we first discuss the twist-3 contributions in the PQCD 
approach. As an illustration, we take the integrand of twist-2, the hard gluon exchange in 
the B meson side, to be 

jtw2 = ^(x 2 )M 2 

' \x x x 2 M 2 B + \k 2X - k 1L \ 2 ][x 2 M 2 + \k 2± \ 2 ] 
where the first term in the denominator is the propagator of the exchanged hard gluon while 



the second one is that of the internal b-quark. As studied in Ref . [[21]] , introducing k± degrees 
of freedom will bring large double logarithms of a s In 2 (k±/ Mb) through radiative corrections. 
In order to improve the perturbative expansion, these effects should be resummed, called k± 
resummation [Ml R2[. Consequently, the Sudakov form factor introduced will suppress the 



region of k\ ~ 0(A 2 ). According to the analysis of Ref. |[34|| , via the Sudakov suppression 



the average < k 2 ± > is of O(AMb) for Mb ~ 5 GeV. Hence, with including k± resummation 
effects, Eq. (Q) becomes 

j tW 2 = ^{X2)M 2 B 

[ Xl x 2 M 2 B + 0(AM B )}[x 2 M B + O(AMb)} ' 1 ' 

Since the fraction momentum x\ is of 0(K/M B ) and 0^-»(x 2 ) oc x 2 {l — x 2 ), it is easy to see 
that at the end-point I tw2 behaves like 

AM B 1 ; 

On the other hand, the integrand of twist-3 is expressed as 

j tW 3_ r K ^ K ,{x 2 )M 2 ^ 



[ Xl x 2 M 2 + 0(AM b )}[x 2 M 2 + 0{AM B )] 



From Ref. [pq] , we find that the twist-3 wave function (j) a K , at the end-point is a constant so 
that 



Ttw3 r K * _ M K * 1 



A 2 A AM B ^ 
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as x 2 — > A/Mb- Hence, the power behavior of I tw3 in Mb is the same as that of I tw2 . We 
note that since the twist-3 one contains the most serious singularity (linear divergence), the 
contribution from a higher twist wave function, such as that of twist-4, should be the same as 
that of twist-3 at most. However, by the definition of the twist wave function, we know that 
the twist-4 one is associated with a factor of r|-*, and its contribution should be one power 
suppressed by than that of twist-3 so that it belongs to a higher power contribution in 
our consideration. In our analysis, its effect will be neglected. 

2.2 Form Factors 



We parametrize the B — > K* transition form factors with various types of interacting vertices 
as follows: 



(K*(p 2 ,e)\V, B( Pl ] 
(K*(p 2 ,e)\A^\B( Pi : 

(K*(p 2 ,e)\T, u q» B( Pl ] 



(K*(p 2 ,e)\T 5 ^ B( Pl 



M B + M K * 



2M K *A {q 2 ) e —^-q, + (M B + M K ») A x {q 2 ) [ ej - ^q, 



~A 2 {q 2 



e -q 



P„ ~ ^J-qp 



M B + M K * V M q 2 



-iT 1 (q 2 )e (ia p p e* a P' 3 q p , 

T 2 (q 2 ) faP -q-e*- qP») + T 3 (q 2 )e* -qL- jf-P,) (10) 



where P = p x + p 2 , q = pi - p 2 , Vp = s 7 M b, = s 7^75 b, T^ u = s ia^b, and 
T^ u = s 20^75 b. According to the PQCD factorization formalism shown in Eq. (|ip, the 
components of form factors defined in Eq. (|1(]) up to twist-3 wave functions are given by 

V{q 2 ) = (1 + r K *)8nC F M 2 / [dx] / b x db x b 2 db 2 4> B {x u h) 

Jo Jo 

x { L£. (x 2 ) - r K * (x 2 (f) v K * (x 2 ) - (- + x 2 )(j) a K * (x 2 )) 

l L fj J 



xE(t ( ^ ) )h(x 1 ,x 2 ,b 1 ,b 2 



+r K * 



(j) v K ,(x 2 ) + r K *^2)]E(t^)h(x 2 ,x 1 ,b 2 ,b 1 )}, 



A (q 2 ) = 8nC F M B [dx] b 1 db 1 b 2 db 2 (f> B {xi, b x ) 

Jo Jo 

x{[(l + rjx 2 )(j) K *(x 2 ) +r K *((l - 2x 2 )(p t K ,(x 2 ) 
+ (- - 1 -2x 2 )(f> s K ,{x 2 ))]E{t^)h{x 1 ,x 2 ,b 1 ,K 

T] 

+2r K *(f) s K 4x 2 )E(t^)h(x 2 , x x , 62, 61)} , 



'.12) 



Mi 2 ) 



lnC F M B ri 



^ [\dx] [°° b 1 db 1 b 2 db 2 <M(xiA) 
Jo Jo 



X 



{[^(x 2 ) +r^((- + x 2 )0^(x 2 ) -x 2 ^.(x 2 )) 

L L V 

xE(t^)h(x 1 ,x 2 ,b 1 ,b 2 ) 

+r K * [<Pl-,(x 2 ) + r K ^2)]E(tfy)h(x 2 ,x 1 ,b 2 ,b 1 )}, 



l+r K 



~-^?>2r K *nC F M 2 B [\dx] [°° b x db x b 2 db 2 <M^iA) 
rj Jo Jo 

x((- - 1)^x2) - (- - l)0^(a: 2 )) J E(t«)/ i (a; 1 ,x 2 ,6 1 ,6 2 ) , 
77 77 



(•1 /"OO 

T\(g 2 ) = 8ttC f M b [dx] b x db x b 2 db 2 fo^iM) 

Jo Jo 

:{[(! + r7x 2 )0^, (x 2 ) + 7*. ((1 - 2x 2 )<f) v K , (x 2 ) 



+(- - 1 - 2x 2 )<j) a K «{x 2 ))} E(tW)h( Xl , x 2 , 61, 62) 
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0^(a: 2 )+#^(x 2 )] J E(t( 2 ))/ i (x 2 ,x 1 ,6 2 ,6 1 )}, 



/■l /"OO 

T 2 (g 2 ) = 8ttC f M1 / [dx] / b 1 db 1 b 2 db 2 (t>B{xiM) 

Jo Jo 

x { [(1 + r]x 2 )4> T K , (x 2 ) + r K * (2 - 77(1 + 2x 2 ))<j> v K , (x 2 ) 
+r K *ri(l -2x 2 ) ( f) a K*(x 2 )}E(tV)h(x 1 ,x 2 ,b l ,b 2 ) 
+r K *V [<Pk* M + (f> a K , (x 2 )] E(t^)h(x 2 , Xl , 62, 61) }, 



rl /*oo 

Tziq 2 ) = 8irC F M B / dx\dx 2 \ bidbib 2 db 2 &i) 

jo Jo 

2 

x { [(1 + r7x 2 )0^« (x 2 ) + r x * (--(* + 2a: 2 ))^. (x 2 ) 

+r x *((l - 2x 2 ))0^(x 2 ) - -^-((I>k*M 
2 

+^(-+^2)^(^2) - ^X*^20^.(^2)) ^^K^l^AA) 



<t> V K .(x 2 )+<i> a K .(x2) - —<p S K *(x 2 )]E(tW)h(x 2 ,X 1 ,b 2 ,b 1 )} 



The evolution factor is given by 



E (t) = a s (t) exp (S B (t) - S K * (t)) , 



where explicit expressions of the Sudakov exponents Sb(k*) can be found in Ref. [Q. The 
hard function of h is written as 

h(x 1 ,x 2 ,b 1 ,b 2 ) = S t (x 2 )K (^/x 1 x 2 i]M B b 1 ) 

x - b 2 )K (y/x^jMeh) I (V^jM B b 3 ) 

+6{b 2 - h)K {^qM B b 2 ) I {y/xwM B bx)] (19) 

where the threshold resummation effect is described by [|34 



91+2CTV3 , c \ 

St(x)= 2 ; [x(l-x)] c . 

V7T1 (1 + Cj 



The hard scales t^ 1 ' 2 ) are chosen to be 



*W = max( v /M 2 r/x 2 , 1/6 1; l/6 2 ) 



t (2) = max( v /M 2 ^ 1 , l/6i, l/6 2 ) . 
For the X* meson distribution amplitudes, we adopt the results given in Ref. |35[ and 



explicitly we have 

4> K *(x) = -^=s(l -x)[l + 0.57(1 -2x) + 0.1(5(1 -2s) 2 -l)], 
V 2^V C 

= ^^[°- 3 ( 1 - 2x )( 3 ( 1 - 2a; ) 2 + 10 ( 1 - 2a; )- 1 ) 

+0.06(1 - 2s) 2 (5(l - 2s) 2 - 3) + 0.21(3 - 30(1 - 2s) 2 + 35(1 - 2s) 4 ) 

+0.36(1 - 2(1 - 2x)(l + ln(l - s)))], 

<^,(s) = - 3 ^[(l-2x)(l + 0.2(1 -2s) + 0.6(10s 2 -10s+l)) 

-0.4s(l - x) + 0.12(1 - 6x - 2 ln(l - a;))], 
3/-T 

</£,(s) = -t=z(1 -x)[l + 0.60(1 -2s) + 0.06(5(1 -2s) 2 - 1)], 

= ^=[^(1 + (l- 2x) 2 + 0.44(1 - 2x) 3 ) + 0.20(3(1 - 2x) 2 -l) 
+0.11(3 - 30(1 - 2s) 2 + 35(1 - 2s) 4 ) + 0.48(2s + ln(l - x))), 
(j) a K *(x) = 3 ^L [(1 - 2s)(l + 0.19(1 - 2s) + 0.81(10s 2 - 10s + 1)) 

-1.14s(l - s) + 0.16(1 - 6s - 2 ln(l - s))] . (20) 

From Eqs. ([i4|)-(|i~6|), at q 2 = we obtain the identities 

A 2 {0) = (l + r K .) 2 ^i(0)-2r^(l + r^)^o(0), 

7\(0) = T 2 (0), (21) 



which are consistent with the leading order model-independent relation T^, |3^, |37|. [IB]. [ID], [10 



A 2 (0) 



l + r K * 



A 1 (0) 



2r 



A'* 



-Mo) ■ 



(22) 



1 - r K * 1 - r K * 

We note that due to the parametrization of Eq. (pi]), there are terms proportional to r 2 K , 



in Eqs. (|TT|) , flT3| ) and (|14l) . In order to guarantee that only tk* dependence appears in the 
left-handed sides of Eq. (|10D, those with r\* should not be dropped. 

3 Angular distributions and physical observables 
3.1 Effective Hamiltonians and Decay Amplitudes 

The effective Hamiltonians of b — ► s £ + £~ are given by pi 



n 



GpctXt 



71 



2^ 



(23) 



with 



H 



H2/j, 

If 

L 5 M 



C 9 (/i) S7 A4 (/i)Pjr,6 

CwSj^PLb , 
If I, 



2rriv 



C 7 (n)sia^q u P R b. 



(24) 



where A t = V t bV t * and Cg(/z), C w and C 7 (/i) are the Wilson coefficients (WCs) and their 
expressions can be found in Ref. for the SM. Since the operator associated with Ciq is 
not renormalized under QCD, it is the only one with the /i scale free. Besides the short- 
distance (SD) contributions, the main effect on the branching ratios comes from cc resonant 
states such as , etc., i.e., the long-distance (LD) contributions. In the literature [12 



H3L [B], |35], fl6j, it has been suggested by combining the FA and the vector meson dominance 
(VMD) approximation to estimate LD effects for the B decays. With including the resonant 
effect (RE) and absorbing it to the related WC, we obtain the effective WC of Cg as 



C e J f 



C 9 (fi) + (3C 1 (fi) + C 2 (fi))\h(x,s) + — J2 k 



7iT(j -> l + l~)M j 



(25) 



^ ' J a 2 — M 2 + ?'MT • 
where h(x, s) describes the one-loop matrix elements of operators 0\ = So-T^-Pl^s c^j^Plc, 



and 2 = s^i^PjJ) c-y^Pic Mj (Tj) are the masses (widths) of intermediate states, and 
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the factors kj are phenomenological parameters for compensating the approximations of 
FA and VMD and reproducing the correct branching ratios Br (B — > J/^/X — > l + FX) = 
Br (B — > J/tyX) x Br (J — > For simplicity, we neglect the small WCs and take 

kj = — 1/ (3Ci (/i) + C*2 (a*))- It is clear that the uncertainty related to this assumption can be 
large [§7| . Moreover, it is questionable whether one can include both quark- level calculations 
with cc-loop and resonances in Eq. Q25p. However, since we are only interested in physics 
behind the various observables at the large recoil we shall not discuss the uncertainties arising 
from Eq. (p5|) . 

Combining Eqs. ([II]) and fl23p, the transition amplitudes for B — > K*£ + £~ (£ = //, e) 
can be written as 

A4*! = f^«>L» + <>i>} (26) 

with 



= iM AH /a/se w (A)PV + 'i2e*(A) + /i3e*-gP / 
= fciW(T(A)PV + 0i<; 



= ijiW^W^V + ^+flse*-?^, (27) 



l^ 9 (g 2 ) 2m b 

2w/, 



to = ' + ^^)Tl{q 2 ) 

mB + rax* q 



h 2 = _( mB + m ^)^(g 2 )-^P-gT 2 V), 

m^e + m^:* q z K F ■ q 

a - C V{q2) 

91 — wo ; ) 

m B + m K * 

92 = -C w (m B + m K *)Ai(q 2 ) , 

<?3 = C W -^-, (28) 
mg + rriK* 

where the form factors associated with superscripts denote the relevant WCs convoluted 
with hard amplitudes and wave functions, described by 



F\q 2 ) = j[dx}[dbW K *{x 2 M)C k {t)T,{{x},{b},M B )^ (29) 
It is worth to mention that for the convenient in the PQCD formalism, Eq. (|29|) can be 



written as F k m Cfc(io) F with F being the form factor and to = \/AMb. 



9 



3.2 Angular Distributions 

In the literature, there are a lot of discussions on B — > K*£ + £~ decays. However, most 
of them have been concentrated on the differential decay rates and lepton polarization and 
forward-backward asymmetries. It is known that the differential decay rates have large 
uncertainties from not only hadronic matrix elements but also the parametrizations of LD 
effects, and the lepton polarization asymmetries are hard to be observed due to the difficul- 
ties of measuring lepton polarizations. Therefore, to test the SM and search for new physics 
in B — > K*£ + £~, it is necessary to find some other physical observables which have less 
theoretical uncertainties but measurable experimentally, similar to the zero positions in the 



forward-backward asymmetries [|38| , |48fl . It is found that if one considers the decay chain of 
B — > K*£ + £~ — > (Kir)£ + £~ , via the study of different angular distributing components, we 
can analyze (a) contributions from both longitudinal and transverse parts of the K* polar- 
ization, (b) T-even (CP conserved) effects from the mixings of longitudinal and transverse 
polarizations of K*, and (c) T-odd effects from the mixings in (b). We note that some 
T-odd effects are suppressed in the SM and thus, measuring these effects could indicate CP 
violation from new physics ]49| . 



To understand dynamical dependence in T-odd terms of e lLVa fi ( f'^ fV Wp? P 13 ■> ft is i n_ 
evitable to investigate the processes of B — > K*£ + £~ — > (Kti)£ + £~ so that the polariza- 
tion A and A' in the differential decay rates, written as dT oc H (A, A') Ai^l Ai^* with 
H(X,X') = e(A) ■ px e*(A') ■ Pk, can be different. From Eq. (EH), we see that M.^} only 



2/U 

depends on C\q. Clearly, T violating effects can not be generated from A^-Mw > but in- 
duced from Mi^MyJ^ as 

wellasA^A^? 1 . Thi s 

can be understood as follows: firstly, for 
the M^M^TrL^L"' part with TrUU' ~ {p^+ +Pe-Pe+ - 9^' Pi- -Pe+), the relevant 
T-odd terms can be roughly expressed by 

M^MffiTrlflS' oc Z 1 Imh 1 h* 3 e(0)-qe^ a ^e w {±)pf + P^ 

+Z 2 Imh 1 h* 2 e(0) ■ p t+ e^ aP q»e* v {±)p a i+ pP 
+Z 3 Imh 1 h* 2 e(T) ■ Pi + e flua/3 q^* u (±)p^ + P 13 (30) 

where Zi (i = 1, 2, 3) are functions of kinematic variables. From Eq. (0), one gets lmh\h\ ~ 
Imhih* 3 ~ ImCg^ (n)Cr(n) . We note that, as shown in Eq. (|3~0D, the T-odd observables 
could be non-zero if the processes involve strong phases or absorptive parts even without CP 
violating phases. By means of Eq. (|25|), Cg^(^) includes the absorptive parts such that the 
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results of Eq. (0) do not vanish in the SM. Secondly, for M^M^TrL^L^' , one gets 



{M { $M%} ] + M^M%) ] )TrL»L^' cc (Imh 2 g* 3 - Imh 3 g* 2 )e, ua pq^(±)p? + pf 3 (31) 



where TrL^L 5 ^' = —Aie^ 1 ' a/5 'q a pi+ p has been used. From Eq. (p8f), we find that Imh 2 g 3 — 
Imh^g 2 is only related to ImCi{ii)Cl and the dependence of ImCg(fi)Ci is canceled in 
Eq. (|3"T|). For the decays of b — > s£ + £~~, since the absorptive parts in Cj(fi) and C±o are not 
expected, a non- vanishing value of ImCi(fi)C^ indicates pure weak CP violating effects. 

In order to derive the whole differential decay rates with the K* polarization, we choose 
that K* helicities are e (0) = (\p K *\,0,0,E K .)/M K * ande(±) = (0, 1, ±i, 0)/\/2, the positron 
lepton momentum p e + = sin^, 0, cos^)/2 with E K * = (Mf - M K * - q 2 )/2y^q 2 



and \pk* \ = y E K , — M K , in the q 2 rest frame, and the K meson momentum pk = 

(1, sin 9k cos0, sin 6^ sin 0, cos9k)Mk* /2 in the K* rest frame where <fi denotes the rela- 
tive angle of the decaying plane between Kit and £ + £~. From Eq. (^), the differential 




decay rates of B — >• K*£ + £ 
dT 



(K7r)£ + £ as functions of angles 9g, 9k and are given by 



d cos 0Kd cos 9id(j)dq 2 




i=l,2 



+ sin 2 0^(1 + cos 2 e t ) (\Mt\ 2 + \Mr\ 2 ) 



-sin 29 K sin 20 t cos0 £ #e(-K + + MT)M° 



+ sin0 J] /m(.M+ - 



-2 sin 2 0^ sin 2 0^ cos 20 ^ #e(.M+.Mr*) + sin 20 ^ Im{MfMi*) 




-2 sin 2 0* cos0^2i?eA^A^* -2i?e.A/^.M 2 * 

+2 sin 20^ sm9t[co8(t){ReMl{Mt* - M 2 *) + Re(Mt - M^)M° 2 ) 

+ sin <P(ImM\(Mt* + -M 2 *) - Im(Mt + M^)M° 2 *)] } (32) 



with 



\P\ 



\]E' 2 - M 2 K 




E 



2M B 




(33) 
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where a = 1(2) while fi = hi (gi) (i = 1, 2, 3). The polarization components M.° a and 
in Eq. (|33|) clearly represent the longitudinal and transverse polarizations, and can 
be easily obtained from Eq. (0), respectively. We note that other distributions for the 
K* polarization and CP asymmetries are discussed in Refs. [5(| and [51 and the photon 
polarization in B — > K*^ — > (Ktt)(1 + 1~) is studied in Ref. [p2|l . 



From Eqs. (|3D|) and (^JJ) , we know that Im{Mt - M~)M°* and Im(M} M~*) are 
from M^M^TrL^L* 1 ', while ImM^Mp + M 2 *) - Im(Mf + Mi)M° 2 * is induced by 
M^Jvi^J^TrL^L?^' . Integrating the angular dependence in Eq. (|32|) , we obtain 

dq 2 1 J 3 x 2% 5 M 2 

x[E(IM°r + lM + i 2 + lM"i 2 )] 

i=l,2 

dT(B K*t 
dq 2 

which conforms the well known equality of 



-Br(K* -> ifvr) 



Br (B -> fT£ + r -»■ Kir£ + r) = Br (5 -»■ K* £ + £~) 5r(iT -»■ Ktt). 



It is interesting to note that by integrating out 6^ and in Eq. (p2[) , we have that 



^r(iT -> K7r){2cos 2 ^ ^ |.M°| 2 



dq 2 dcos6 K 2 10 n 5 m 2 B r _ ] ,, 

+ sin 2 ^ £ {\Mt\ 2 + \MT\ 2 )} , (34) 

i=l,2 

which allow us to define normalized longitudinal and transverse polarizations of K* by 

E i=1 , 2 |-M°l 2 



V L (q 



2\ 



AI2 ' 



= E,E,I^P ' (35) 



respectively. 

3.3 Physical Observables 



From Eq. (|32"D, it is clear that there are 9 different helicity combinations in the amplitudes. 
As we will show next, among them, 6 are T-even and 3 T-odd. If each component can be 
extracted from the angular distribution, we should have 9 physical observables, which can 
be measured separately in B — > K*£ + £~ decays. To archive the purpose, we will propose 

12 



some proper momentum correlation operators, so that each component of Eq. fl3"2|) can be 



singled out and measurable experimentally. In the following discussions, we use the K* rest 
frame. The coordinates of relevant momenta are choosing as follows: 

p B = (7-Mb, 0, 0, -j(3M B ), p £+ = E i+ (p° i+ , sin 6 t , 0,p 3 i+ ), 
p° e+ = 7 2 (1 + (3 2 - 2(3 cos 9t), p 3 £+ = 7 2 ((1 + (3 2 ) cos9 £ - 2(3), 

* = !.7 1 



t? MK * TP /ofiN 

= U K = — — , = W£+ = — - , (36) 

where (3 and 7 are the usual Lorentz transformation factors. From Eq. ([36]), we have 

Pi+ x p B 



sin Qi 
sin 

COS(/> 



Pb x 

(Pb x p K ) ■ (p e + x p B ) 
\Pb x x p B \ 



sm0 = p B — — — . (37) 

\Pb x Pa-||pb x Pi+\ 

To relate the above angles to those in Eq. (p7|), we use momentum correlations denoted by 
Oi and define the physical observables by 

dT 

dq 2 



(d) = I O^iue^v^— (38) 



where Ui(ug K ,vg e ) are sign functions, i.e., 

uJi{ue K ,v ee ) = 1 (39) 

\ue K vg e \ 

with uq. being sin^j or cos^j. The asymmetries Ai and statistical significances £j of Oj are 
given by 

A(a 2 ) = <<Dl> 

£j ( g 2) = <<%> (4Q) 



We can also define the integrated asymmetries and statistical significances by 

J O i uj i (u eK ,ve t )dT 



A- 

~ JdT 

J OiUi{ue K ,v Bl )dT 



JdT-JOfdT 
13 



The numbers of B mesons required to observe the effects at the n a level are given by 



rr 



Ni = o, (42) 

BR- AV 



where Ai represent the asymmetries or the statistical significances and BR the branching 
ratios of B -> Kn£ + £-. 

To study the various parts of the angular distributions in Eq. (|32]) , we use nine operators 
Oi (i — 1, 2, • • • , 9) and sign functions u>i as follows: 



O 



\P£+xPb\ 2 JpbXPk 12 



' \pb\ 2 u 2 + \Pb\ 2 ^k 

lo\ = ui(sm 9k, sin 9i) , (43) 



Q = 2 \Pb x p K \ 2 _ \p e + *Pb\ 2 
~ 1 \Pb\ 2 u\ \pb\ 2 ^ 2 + 

u 2 = Ui, (44) 



(Pb x Pk) ■ {Pt+ x p B ) 



Wb x p K \\Pi+ x p B \ 

uj 3 = uj 3 (cos 9 K , cos 9i) , (45) 



, [(Pb x pW) • (p i+ x p B )] 2 _ 2 (p B ■ p> x p K ) 2 



Ipb x pV| 2 |pV x p B | 2 |p B x pV| 2 |p£+ x p B \ 2 

uJi = lui , (46) 



o 5 = 1, 

u; 5 = CL> 5 (sin0tf,cos0£) , (47) 
Uq = ujq(cos 9k, sin 6^) , (48) 



n i - i Pk-{pb>< Pi+) 

Oj = \Pb 



\Pb x Pk\\pb x 

c^7 = uj-? (cos 9k, cos 9 e) , (49) 



I _ I (pB Pe+ x p K ) (pb x pW) ■ ($+ x p B ) 
= \Pb\ 



Wb x pV| 2 |p> x p B \ 2 
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UJ 8 



(50) 



o 9 



o 7 

UJ 6 . 



(51) 



It is clear that the first six operators Oi (i = 1 — 6) in Eqs. (f43|)-(f48|) are T-even observables, 
whereas the last three Oj (i = 7 — 9) T-odd ones. We remark that the operators and sign 



functions in Eqs. ([43|)-([51]) are the simplest ones to discuss the momentum correlations. 
From Eqs. ©, ©, and ©-(H), we find that 



A 2 (q 2 

A 4 (q 2 

A 5 (q 2 
A 6 {q 2 
Mq 2 

As(q 2 

A 9 (q 2 



9r 



i=l,2 



64tt 



9r l=h2 

167T 



Ed-^T + I-^ 



9r 



E i?e(>l+ + Mr)^ 



0* 
i ) 



i=l,2 



-2~ E MMtMr) , 

91 i= l,2 

-2 • ^{2ReMtMp - 2ReM^M 2 *) , 
2tt 2 

2 • (# e .M?(.M+* - M 2 *) + ife(A<+ - A^r)^"*) 

61 o 



16/T 



9r 



E Im{Mt - Mj)M 



0* 

i i 



i=l,2 
7T 



91 i= l,2 

2vr 2 

2 ■ (j m M\{Mt + M^*) - Im{Mt + M^Ml*) , 

6L o 



(52) 



where r = Q^/9J2\J2i \Mf \ . We note that the asymmetry .4.1(2) hi Eq. (|52|) is related 
to the longitudinal (transverse) polarization of K* in Eq. Q55|). We can also evaluate £i(q 2 ) 
similar to those in Eq. fl52|) except the denominators due to / OfdT/dq 2 . 



4 Numerical Analysis 

In our numerical analysis, we use /# = 0.19, fx* = 0.21, = 0.17, M# = 5.28, M^* = 
0.892, nib = 4.8 GeV, Xt = 0.04, a em = 1/129, and c = 0.4, and we take the B meson wave 
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Table 1: 


Form factors for B — » 


K* at g 2 = 


in various QCD mo 


dels. 


Model 




V(0) 


4,(0) 


4(0) 


4(0) 


71(0) 


Ts(0) 


LEET | 




0.36 ±0.04 




0.27 ±0.03 




0.31 ±0.02 




QM H 


I 


0.44 


0.45 


0.36 


0.32 


0.39 


0.27 


LCSR [ 


m 


0.399 


0.412 


0.294 


0.246 


0.334 


0.234 


LFQM 


© 


0.35 


0.32 


0.26 


0.23 


0.32 


0.21 


PQCD (I) 
(II) 


0.355 
0.332 


0.407 
0.381 


0.266 
0.248 


0.202 
0.189 


0.315 
0.294 


0.207 
0.193 



function as 

M*> b) = N B x\l - x) 2 exp [ - \C—? ~ 4^1 > ( 53 ) 
where ujb is the shape parameter [53| and Nb is determined by the normalization of the B 



meson wave function, given by 



1 h 



dx(J)B(x,0) = J " (54) 



Since the PQCD can be only applied to the outgoing particle of carrying a large energy, 
where a small coupling constant a s expansion is reliable, we only perform our numerical 
analysis in q 2 < 10 GeV 2 . 

4.1 form factors 

In Table [I], we show the form factors parametrized in Eq. ( |i~0D with (I) uob = 0.40 GeV 
and (II) ujb = 0.42 GeV at q 2 = 0. As comparisons, in the table we also give results from 
the light cone sum rule (LCSR) f38f , quark model (QM) |39|1 , and light front quark model 
(LFQM) pD| . Since in the large energy effective theory (LEET) seven independent form 
factors in Eq. ( |i~0D can be reduced to two in the small q 2 region |36[], in Table [I], we only 
show 7\(0), V(0), and 4(0) [0 extracted by combining the LEET and the experimental 
data on5^ K*j. 

In our following numerical analysis, we only take the minimal results in the LCSR, which 
are consistent with those from the extraction of the LEET, as the representation of the 
LCSR. From Table [3], we find that our results from the PQCD agree with those from all 
other models except the QM. 

It is interesting to point out that the decay branching ratio (BR) of B — > <pK°* is found 
to be 1.7 (1.5) x 10~ 5 for uo B = 0.40 (0.42) |>§, comparing with the recent BELLE's result 
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of 1.3 

i 0-21 x 10 5 [^]. Here, the overwhelming contributions to BRs of B — > <fiK* are 
from the longitudinal parts, where the form factor Aq plays an essential role. We remark that 
A does not appear in our present analysis due to the light lepton mass neglected. However, 
one can obtain the value of A if more accurate measurements on the modes of B — > <j)K* 
are available in the near future. After getting A (0) and ^i(O) , we can find A 2 (0) from 
the identity in Eq. (pij). Furthermore, by using the relations among the form factors in 
the HQET |J7|], one can easily get 23(0) as well. In sum, in terms of the measurements of 
B -> <frK* and B -> K*^ together with the HQET and LEET, all form factors at q 2 = for 
B — > K* can be extracted model-independently. 

In Figures [L]-0, we display the form factors V(q 2 ), A 0t i^(q 2 ), and Ti j2i3 (g 2 ) as functions 
of s = q 2 /M 2 , in the LCSR, QM, LFQM, PQCD (I), and PQCD (II), representing by the 
solid, dash-dotted, dotted, square and circle curves, respectively. 

4.2 Differential decay rates 

We now present the dilepton invariant mass distributions for B —>■ K*l + l~ by integrating 
all angular dependence in the phase space. The distribution for the decay mode with a 
muon pair in various QCD approaches is shown in Figure || where (a) and (b) represent the 
results with and without resonant effects, respectively. From the figures, we find that the 
PQCD results are consistent with the minimal ones in the LCSR approach due to similar 
form factors in the lower q 2 region in both models. Here, we have set that all WCs are 
involved at the m& scale for all QCD approaches except the PQCD one. As emphasized in 
Sec. Ill, in the PQCD formalism WCs should be convoluted with hard parts and meson 
wave functions of B and K*. Due to the hard gluon exchange, with the momentum squared, 
|/cjj 2 , being off-shellness in magnitude of 0(AMg), dominated in the fast recoil region, the 
PQCD approach involves a lower scale |2B], [EJ. As a consequence, even using the concept of 
the naive factorization, where the decay amplitude is expressed by the product of the WC 
and the corresponding form factor, the typical scale to m the PQCD should be much less 
than Mb or Mb/ '2. To illustrate the scale dependent on the WCs, we display the Cr(ji) and 
Cg(/f), renormalized by themselves at m&, as functions of //-scale in Figure In Figure [TU], 
we show the decay rate of B — > K*fi + fi~ with the relevant WCs fixed at to — 1-3, 1.5, 1.7 
and 5.0 GeV, respectively. From the figure, we find that the result with t = 1.3 GeV is 
compatible with that from the formal PQCD approach. Similar conclusion is also expected 



17 



for the electron mode. 



4.3 Physical observables 

Because the numerical values of Ai(q 2 ) are similar to Ei(q 2 ), in the following numerical 
calculations, without loss of generality we concentrate on Si{q 2 ). Moreover, we will not 
discuss the contributions from ©3,4,7,8 since they are very small. 



In Figures [XT] and [T2|, we show the statistical significances of 0\ and O2, related to the 
longitudinal and transverse polarizations of K* , in various QCD approaches, respectively. 
From these figures, we see that the differences among different QCD approaches are in- 
significant, i.e., they are not sensitive to hadronic effects so that they can be used as good 
candidates to test the SM as well as search for new physics. 

In Figures [13] and |14], we display £5(3) and Eq(s) as functions of s — q 2 jM\, which 
correspond to the angular distributions of cosde and sin 29k in Eq. ( |32|) and depend both 
on ReCg^C* and ReC 7 C± , respectively. From the figures, we see that the zero points of e 5 
and Eq in the PQCD are quite different from others. Hence, by measuring these distributions, 
especially those zero points, we can distinguish the PQCD results from other QCD models. 

To show T violating effects, we concentrate on the T-odd operator of Og and consider new 
CP violating sources beyond the CKM. In the SM, the contribution to eg is less than 0(1%). 



As illustrations, in Figures [15] and |16|, we present our results by taking (i) ImCj = 0.25 
and (ii) ImC-j = 0.25 and ImC\Q = —0.20 with the others being the same as those in 
the SM, respectively. One possible origin of having these imaginary parts is from SUSY 
where there are many CP violating sources. It is interesting to see that the CP violating 
effect in Figure [16| can be as large as 30% in these models with new physics. We emphasize 
that a measurement of such effect is a clear indication of new physics as contrast with the 
decay rates for which one could not distinguish the non-standard effect due to the large 
uncertainties in various QCD models as shown in Figure |8]. Finally, we note that unlike eg, 
£7,8 receive contributions from the absorptive parts in Cg^(n) in the SM and they conserve 
CP. On the other hand, they are much smaller than eg in new physics models such as the 
ones in (i) and (ii). Due to the uncertainty in the form of Cg^(fi) in Eq. ([25]), we shall not 
discuss them further. 
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5 Conclusions 



We have studied the exclusive decays of B — > K*£ + £~ within the framework of the PQCD. 
We have obtained the form factors for the B —>■ K* transition in the large recoil region, 
where the PQCD for heavy B meson decays is reliable. We have found that the form factors 
at q 2 = are consistent with those from most of the other QCD models, in particular, the 
LEET combined with the HQET and the experimental data on B — > K*j. We have related 
the angle distributions in the decay chains of B — > K*(K7r)£ + £~ with 9 physical observables 
due to the different helicity combinations in B — > K*£ + £~ . In particular, we have shown 
that the T-odd observable can be used to test the SM and search for new physics, which is 
< 0(1%) and up to 0(10%), respectively. Finally, we remark that to measure such an 10% 
CP violating effect experimentally at the 1 a level, for example, in B — > Kirfi + jjT one need 
at least 7.7 x 10 7 BB, which is accessible in the current B factories at KEK and SLAC. 
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Figure 1: Form factor for V(q 2 ) as a fuction of g 2 /Mj>. The curve with squares (circles) stands 
for the PQCD calculation with ujb = 0.40 (0.42) and the solid curve denotes the minimal 
values in the LCSR [|38|, while the dash-dotted and dotted ones represent the results of the 



QM and LFQM, respectively. 
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Figure 7: Same as Figure [l] but for T 3 (q 2 ). 




Figure 8: Differential decay rate of B — > K*fi + fi~ as a function of s = q 2 : /M|. The solid, 
dotted, and dashed curves stand for the results of the LCSR, QM, and LFQM; the upper 
and lower dash-dotted ones are those from the PQCD (I) and (II), and (a) and (b) represent 
the results with and without resonant effects, respectively. 
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0.75 




Figure 9: Wilson coefficients as a function of \x normalized by themselves at fi = m&. The 
solid and dashed curves are for a 7 = C 7 (fi)/C 7 (m b ) and a 9 = C 9 (/z) / C 9 (mb) , respectively. 




Figure 10: Differential decay rate of B — > K*[i + fi~ for different \x scales. The dotted, 
dashed, dash-dotted, double-dot-dashed curves stand for fi = 1.3, 1.5, 1.7, and 5.0 GeV, 
respectively, while the solid one expresses the result of the full PQCD formalism. 
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Figure 11: Statistical significance Si(q 2 ) of 0\ as a function of s = q 2 /Mg. The solid, dotted, 
dashed, and dash-dotted curves stand for the results of the LCSR, QM, LFQM, and PQCD 
(I), respectively. 
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